High-order extension of the recursive regularized lattice Boltzmann method

This thesis is dedicated to the derivation and the validation of a new collision model as a stabilization technique for high-order lattice Boltzmann methods (LBM). More specifically, it intends to stabilize simulations of: (1) isothermal and weakly compressible flows at high Reynolds numbers, and (2) fully compressible flows including discontinuities such as shock waves. The new collision model relies on an enhanced regularization step. The latter includes a recursive computation of nonequilibrium Hermite polynomial coefficients. These recursive formulas directly derive from the Chapman-Enskog expansion, and allow to properly filter out second- (and higher-) order nonhydrodynamic contributions in underresolved conditions. This approach is even more interesting since it is compatible with a very large number of velocity sets. This high-order LBM is first validated in the isothermal case, and for high-Reynolds number flows. The coupling with a shock-capturing technique allows to further extend its validity domain to the simulation of fully compressible flows including shockwaves. The present work ends with the linear stability analysis(LSA) of the new approach, in the isothermal case. This leads to a proper quantification of the impact induced by each discretization (velocity and numerical) on the spectral properties of the related set of equations. The LSA of the recursive regularized LBM finally confirms the drastic stability gain obtained with this new approach

Extension d'ordre élevée pour les méthodes Boltzmann sur réseau régularisées par récurrence

Ce manuscrit est consacré au développement et à la validation d'un nouveau modèle de collision destiné à améliorer la stabilité des modèles lattice Boltzmann (LB) d'ordre élevés lors de la simulation d'écoulements : (1) isothermes et faiblement compressibles à nombre de Reynolds élevés, ou (2) compressibles et comprenant des discontinuités telles que des ondes de choc. Ce modèle de collision s'appuie sur une étape de régularisation améliorée. Cette dernière inclut désormais un calcul par récurrence des coefficients hors-équilibre du développement en polynômes d'Hermite. Ces formules de récurrence sont directement issues du développement de Chapman-Enskog, et permettent de correctement filtrer les contributions non-hydrodynamiques émergeant lors de l'utilisation de maillages sous-résolus. Cette approche est d'autant plus intéressante quelle est compatible avec un grand nombre de réseaux de vitesses discrètes. Ce modèle LB d'ordre élevé est validé tout d'abord pour des écoulements isothermes à nombre de Reynolds élevé. Un couplage avec une technique de capture de choc permet ensuite d'étendre son domaine de validité aux écoulements compressibles incluant des ondes de choc. Ce travail se conclut avec une étude de stabilité linéaire des modèles considérés, le tout dans le cas d'écoulements isothermes. Ceci permet de quantifier de manière distincte l'impact des discrétisations en vitesse et numérique, sur le comportement spéctrale du jeu d'équations associé. Cette étude permet au final de confirmer le gain en stabilité induit par le nouveau modèle de collision.This thesis is dedicated to the derivation and the validation of a new collision model as a stabilization technique for high-order lattice Boltzmann methods (LBM). More specifically, it intends to stabilize simulations of: (1) isothermal and weakly compressible flows at high Reynolds numbers, and (2) fully compressible flows including discontinuities such as shock waves. The new collision model relies on an enhanced regularization step. The latter includes a recursive computation of nonequilibrium Hermite polynomial coefficients. These recursive formulas directly derive from the Chapman-Enskog expansion, and allow to properly filter out second- (and higher-) order nonhydrodynamic contributions in underresolved conditions. This approach is even more interesting since it is compatible with a very large number of velocity sets. This high-order LBM is first validated in the isothermal case, and for high-Reynolds number flows. The coupling with a shock-capturing technique allows to further extend its validity domain to the simulation of fully compressible flows including shockwaves. The present work ends with the linear stability analysis(LSA) of the new approach, in the isothermal case. This leads to a proper quantification of the impact induced by each discretization (velocity and numerical) on the spectral properties of the related set of equations. The LSA of the recursive regularized LBM finally confirms the drastic stability gain obtained with this new approach

Extensive analysis of the lattice Boltzmann method on shifted stencils

Standard lattice Boltzmann methods (LBMs) are based on a symmetric discretization of the phase space, which amounts to study the evolution of particle distribution functions (PDFs) in a reference frame at rest. This choice induces a number of limitations when the simulated flow speed gets closer to the sound speed, such as velocity-dependent transport coefficients. The latter issue is usually referred to as a Galilean invariance defect. To restore the Galilean invariance of LBMs, it was proposed to study the evolution of PDFs in a comoving reference frame by relying on asymmetric shifted lattices [N. Frapolli, S. S. Chikatamarla, and I. V. Karlin, Phys. Rev. Lett. 117, 010604 (2016)]. From the numerical viewpoint, this corresponds to overcoming the rather restrictive Courant-Friedrichs-Lewy conditions on standard LBMs and modeling compressible flows while keeping memory consumption and processing costs to a minimum (therefore using the standard first-neighbor stencils). In the present work systematic physical error evaluations and stability analyses are conducted for different discrete equilibrium distribution functions (EDFs) and collision models. Thanks to them, it is possible to (1) better understand the effect of this solution on both physics and stability, (2) assess its viability as a way to extend the validity range of LBMs, and (3) quantify the importance of the reference state as compared to other parameters such as the equilibrium state and equilibration path. The results clearly show that, in theory, the concept of shifted lattices allows the scheme to deal with arbitrarily high values of the nondimensional velocity. Furthermore, just like the zero-Mach flow for the standard stencils, it is observed that setting the shift velocity to the fluid velocity results in optimal physical and numerical properties. In addition, a detailed analysis of the obtained results shows that the properties of different collision models and EDFs remain unchanged under the shift of stencil. In other words, by introducing a velocity shift in the stencil, the optimal operating point, in terms of physics and numerics, will also be shifted by the same vector regardless of the EDF or collision model considered. Eventually, while limited to the D2Q9 stencil with the nine possible first-neighbor shifts, the present study and corresponding conclusions can be extended to other stencils and velocity shifts in a straightforward manner